ANALYTICAL CALCULATION OF THE NEUTRONS SPECTRUM 
FOR DIRECT MEASUREMENT OF N-N SCATTERING 
AT PULSED REACTOR YAGUAR 
V. K. IgnatovictQ 
FLNP JINR Dubna 141980 RF 

Abstract 

Analytical calculation of a single neutron detector counts per YAGUAR reactor pulse 
is presented and comparison with coincidence scheme is given. 
PACS: 28.20.-v; 28.20.Cz; 29.30.Hs; 14.20.Dh; 25.40.Dn 

1 Introduction 

There is a project to measure directly n-n collision for checking charge symmetry of nuclear 
forces pQ. It is accepted that the best neutron source to perform such measurements is the 
Russian pulsed YAGUAR reactor. Some preliminary measurements and numerical simulations 
for expected experimental geometry had been performed [2]. We want to show here an analytical 
approach to calculations. First we obtain analytical momentum spectrum of scattered neutrons, 
then the time of flight spectrum of neutrons detected by a single counter. After that we consider 
coincidence scheme where we have two detectors, and calculate time of flight spectrum for one 
detector and delay time spectrum for the second one. We considered coincidence scheme because 
from the very beginning of discussions about the project, and all the time during preparation of 
the experiment, many people continue to express the opinion that the coincidence scheme has 
an advantage comparing to the single detector measurement. They claim that loss of intensity, 
which they usually estimated at the level of 20%, will be surpassed by much higher suppression 
of background. We show here analytically that in the coincidence scheme effect is so much 
suppressed, that the question about the background level becomes irrelevant. 

2 Estimation of the effect 

The scheme of the experiment is presented in Fig. 1 borrowed from [1] . The YAGUAR reactor 
1 gives a pulse of length t p = 0.68 ms, during which a huge amount of neutrons with flux density 
$ = 0.77 x 10 18 n/cm 2 s is released. After a moderator at room temperature T neutrons in the 
thermal Maxwellian spectrum arrive at the volume 2 (V = 1.13 cm 3 ), where they collide with 
each other and some of them after collision fly along the neutron guide 3 with collimators 4, 
and arrive at the detector 5, where they are registered with ~ 100% efficiency. The collimators 
4 determine the solid angle AQ = 0.64 x 10~ 4 , at which the volume V is visible by the detector. 
The estimated number of neutrons that can be registered at a single pulse is equal to 

N e = 2n 2 Vt p v T \b\ 2 dVt ) (1) 

where factor 2 takes into account that the detector can register scattered neutron or neutron- 
scatterer. The square of the scattering amplitude \b\ 2 is defined as: \b\ 2 = |6 | 2 /4, where bo is the 
singlet scattering amplitude, which is accepted to be 18 fm, and factor 1/4 is statistical weight of 
the singlet scattering. Therefore \b\ 2 = 8.1x 10~ 25 cm 2 . The speed vt corresponds to the thermal 
speed Vt = 2200 m/s, and the factor fr|^| 2 determines number of collisions in the neutron gas 
per unit time. The factor n 2 is the square of the neutron density: n = $/t>T = 3 x 10 12 cm~ 3 . 
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Figure 1: Scheme of the experiment on direct measurement of n-n scattering pQ. 1 — reactor 
core; 2 — volume of collisions; 3 — neutron guide; 4 — collimators; 5 — detector; 6 — neutrons 
trap. 

After substitution of all the parameters into ([I]) we find N e pa 170 neutrons per pulse. However 
it is the estimation number. To find real number counted by the single detector, N s , it is 
necessary to calculate the scattering process. Calculation shows that N s = FN e , where factor 
F is of the order unity. Monte Carlo calculations in [Tj give F = 0.83. Analytical calculations 
presented below give F = 0.705. The number of neutrons per pulse counted at coincidence, if 
the neutrons trap 6 is replaced by another detector, can be estimated as 

N ec = N s dnr/t T , (2) 

where r is the width of the coincidence window, tj< = L/vt is the average length of measurement 
time after the reactor pulse, and L pa 12 m is the average distance between collision volume 
and the detectors. In the experimental scheme of Fig. 1 the time is of the order 5 ms. If we 
accept t pa t p = 0.5 ms, then the ratio r/tx is 0.1. The factor <El is included in fl2]), because 
only neutrons in this solid angle will be registered by the second detector. The total factor, 
which suppresses the estimated number of neutrons registered per single pulse in coincidence 
scheme, is of the order 10~ 5 , therefore the estimated number of counts in coincidence scheme 
will be 10~ 3 , so the experiment becomes non feasible, and the level of the background, which is 
determined by neutron scattering on the residual gas atoms present at even very good vacuum 
conditions, becomes irrelevant. The analytical calculations, presented below, show that the 
real number of counted neutrons in coincidence scheme contains even additional small factor 
F c = 0.15. 

3 The analytical calculation of neutron scattering in the 
thermal neutron gas 

Our calculations will be based on the standard scattering theory of neutron scattering in the 
atomic gas. Our main feature is that we shall make calculations directly in the laboratory 
reference frame without transition to the center of mass system. First we remind all the 
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definitions of the standard scattering theory and then present analytical calculations of all the 
required integrals. 



3.1 The standard scattering theory 

The standard scattering theory starts with the Fermi golden rule, according to which one can 
write down the probability of the neutron scattering per unit time on an arbitrary system as 

dw(k, t -> k f , Xi -> X f ) = y \{X f , k f \U\\, k t )\ 2 6{E fk + Efx - E ik - E iX )p(E fk ), (3) 

where |fcj > , |A, > are initial, \kf >, |A/ > are final states of the neutron and system with 
energies E ik , E iX , Ef k , Ef X respectively, U is the neutron-system interaction potential, which 
in the neutron atom scattering is accepted in the form of the Fermi pseudo potential 

h 2 

U= — 47rW(ri-r 2 ). (4) 

Here r±, r 2 are positions of the neutron and the system, p(Ef k ) is the density of the neutron 
final states 

p{E k ) = [^)d% (5) 

E k = h 2 k 2 /2m, m is the neutron mass, and L is the size of some arbitrary space cell. 

We suppose that the system is an atom with mass M = m, and momentum p. The initial 
and final states of the neutron and atom are described with similar wave functions 
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/ >= jjj5 e M ik ij r ), Wf>=\Pi,f>^j^ e M i Pij r ), ( 6 ) 



where kij and p i j are initial and final neutron and atom momenta respectively. 
The flux density of the single incident neutron is 

ji = hki/mL 3 . (7) 

The scattering cross section at the given initial and final states is the ratio 

da(ki -> k f , Pi -> p f ) = -dw(ki -> k f , Pl -> p f ). (8) 

Ji 

At the next step we need to sum this cross section over final states of the system and average 
over initial states. In our case summation over the system final states is the integration over 
density of the atomic final states 

P(E pf )={£) d"p f . (9) 

This integration gives the cross section for the given initial states as 

da(ki -> kf,Pi) 

2itm L 9 d 3 k f 
'h 2 ki (2ttY 

where E p = % 2 p 2 /2M, E k = h 2 k 2 /2m, and the additional factor 2 means that the atom and 
neutron are the same particles, therefore we can detect with the same probability the scattered 
neutron in the phase element d 3 kf or an atom in the element d 3 pf. 



= 2- 



— / j d 3 p f (p f , k f \U\ Pi , hi) 5(E fk + E fp - E ik - E ip ), (10) 
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For our experiment we need not a cross section, but the number of the neutrons dN{k il p { , kf) 
scattered in the element d 3 kf. This number is determined by the number of collisions of neu- 
trons with atoms, so the number of scattered neutrons is equal to 

dN(ki,Pi,kf) = dn a (pi)dn n (ki)Vdt p vda(pi, ki — > kf), (11) 

where dn a (pi), dn n (ki) are the number densities of atoms and neutrons with initial momenta 
p i and ki respectively, v — 7i|p i — ki\/m is the relative neutron-atom velocity, and V, dt p are 
volume and time, where collisions create detectable neutrons. 

Since our atoms and neutrons have the same Maxwellian distribution with the temperature 
T, the densities dn a (pi) and dn n {k.i) are 

dn a (q) = dn n (q) = n j^jy/5 exp [ -— I , (12) 

where n is the average neutrons density, the letter T denotes reduced temperature T = 
mfcs[T]/^, 2 , and [T] is the temperature in Kelvin degrees. To find the total number of neutrons 
dN(kj) scattered into element d 3 kf of the final momentum space we must integrate (TTTT) over 
dn a (pi)dn n (ki), after which we get 



xexp(-?iiM) Jd 3 p,\{ P ,,k,\V\p„k i )\'s{k}+p)^kl-pf). (13) 



The matrix element of the potential (HI) is 



T) 2 (2tt\ 3 

(p f , k f \V\ Pi , ki) = Aacb—ijjLsfa + k i -p f - k f ), (14) 



and its square is 



|2 fh 2 \ 2 (2tt) 3 



( Pfl k f \V\ Pt ,h) = |4< I - j ^J-Sfa + h-pf-kf). (15) 

After substitution of ffT5j) into ffl3|) we can extract |6| 2 from the square of the matrix element, 
dQ from d 3 kf and introduce the thermal speed vt = h\r2T/m. As a result we obtain 



dk 

dN(k f ) = N e g{kf)^L, (16) 



where N e is given in ([I]), and <?(&/) is 

« k >> = J f 

x Jd 3 p f e X pL?J 1 ^\s(p i + k,-p f -k l )S(k 2 f +p}-k?-pt). (17) 
Integration over c? 3 pj gives 

*(*/) = / ^/exp (-^S) / d^J^M^fcJ +P ) - k 2 - (P - A,) 2 ), (18) 

where P = pj + k / is the total momentum of two particles. 

With all these definitions in hands we can directly calculate the spectrum of scattered 
neutrons 



4 



3.2 Analytical calculation of the integrals 

First we calculate the integral 

|P-2fei| 



Q(k f , Pf ) = J d% 1 il S(k}+p 2 f - k 2 - (P - k t f) 

= 2 J ^\2k t - P\S((k f - Pf ) 2 + (2k, - Pf). (19) 
After change of variables 2k, — P = u we obtain 

1 r d ii 

Q(kf,P f ) = - 2 j «j^p[*(« a ~ <?l (2°) 

where q 2 = (kf — pf) 2 . 

After representation ud 3 u = (u 2 du 2 / '2) dipd cos 6, where polar axis is chosen along the vector 
P, we can integrate over dtp and d(u 2 ). As a result we get 

nq 2 d cos 9 



1 2W + 2p qa »t + i» - (21) 



-l 



Integration over dcos9 gives 



Q(k flPf ) = ^{q + P-\q-P\). (22) 



The last factor is equal to 2q, if q < P, and it is equal to 2P, if q > P. Which one of these 
inequalities is satisfied depends on the angle Of between vectors kf and Pf. Inequality q < P 
is satisfied, when cos8f > 0, and inequality q > P is satisfied, when cos6f < 0. Therefore Eq. 
(I22jl is representable in the form 

Q(kf,p f ) = ixq (e(cos#/ < 0) + e(cos#/ > 0)j^j , (23) 

where Q(x) is the step function equal to unity, when inequality in its argument is satisfied, and 
to zero in the opposite case. 

3.3 The spectrum of neutrons, counted by a single detector 

Substitution of fl23|) into f|T8|) gives 

g(k f ) = J d 3 pfw(kf,p f ), (24) 

where 

= ^^ ex p {-^rj Q( k f>pf)- ( 25 ) 

To obtain spectrum of neutrons counted by a single detector we represent d 3 pf = p 2 dpfdQf, 
and integrate Q(k,p) over dQf. As a result we obtain (in the following we omit subscripts / 
of variables) 

/■! (h-n\ 2 \ 



I(k,p)= jQ(k,p)dQ = 2n 2 IJ dcos6\k - p\ + J dcosd 



( fc-p) 
\k + p\ 



Figure 2: Spectrum y 2 F(y) exp(— y 2 ) of the neutrons detected by a single detector in dimen- 
sionless units y = v/v T . 



2n 2 (J dcos9\k + p\+J dcos9 



2(k 2 +p 2 ) 
\k +p\ 



\k +p\ 



M 

pk 



(p 2 + k 2 )(p + k-Jp 2 + k 2 ). 



Substitution of fl26l) into (|25|) and change of variables x = p/k, y = k/^/2T gives 



g{k f ) = f(y) = -eM-y 2 )y 2 J(y)-, 



TV 



where 



J( y ) = 2y 4 / 2xdxexp(-x 2 y 2 ){l + x 2 ) \(x + 1) - Vx 2 + 1 
Jo 1 



Integration by parts gives 

roc 

J(y) = 2y 2 / dxex.p(-x 2 y 2 )(l + 2x + 3x 2 - 3xVx 2 + 1) = y0F + My), 
Jo 

where 

J x (y) = 2y 2 / xdxexp(-xV)(2 + 3x - 3VX 2 + 1) 



and 
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-1 + 3 / x<irexp(— x 2 ?/ 2 )(l 
Jo 



^L= ) = _ 1+3 £ { 1-^[1 -.<„)]>, 



2 [V 2 

$(y) = — p= / c?xexp(— x ). 
y 7T Jo 



Substitution of pjj into gives 



7T 



J(y) = ^ - 1 + 3^{1 - eMl - <%)]}• 



(26) 

(27) 
(28) 

(29) 



(30) 
(31) 

(32) 



The momentum spectrum f(y) from Eq. (l2~7l) with account of ( l32l) is shown in Fig. 2. Numerical 
integration of this function gives F = Jq 00 f(y)dy = 0.705. 

3.4 Time of flight spectrum of a single detector 

In the experiment the time of flight (TOF) spectrum is measured. To transform ( J27|) into TOF 
spectrum we multiply it by unity 



1 = dt5(t- L/vry), 



(33) 
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where L is the distance between scattering volume and the detector, and integrate over dy. 
After that we obtain 

L ft \ y 



N s (y(t)) = f(L/v T t)— = f{y) y - (34) 



4 Registration by two detectors in coincidence 

Let's consider the case, when neutrons are registered in coincidence by two detectors on the 
opposite sides of the collision volume. It means that the angle between kf and Pf is approxi- 
mately 180°. Since we register both neutrons, we should not integrate (I2"3"j) over d 3 pf. Instead 
we should accept kfPf < 0, Pf ~ — kf, and d 3 pf = pj-dpfdQ with the same dfl as in d 3 kf. 
Taking into account Eq. (j!6j) . (l24"j) and fl25|) we can represent the number of neutrons counted 
by two detectors as 

2 k^dkf ( p 2 + A; 2 \ 

dN(k f ,p f ) = Njn pf -jJ^qp 2 f dp f exp (^-^^J • (35) 



After transformation to dimensionless variables y = kf/y2T and z = pf/\/2T we get 

dN(k f , p f ) = N e dQ p G{y, z)dydz, (36) 

where 



2^2 



G(y, z ) = 2 -y^{ y + z ) eM-y 2 - ^ 2 ), (37) 



7T Z 

and we replaced q by kf + pf. 

To get TOF spectrum in one detector and coincidence count in the second one with coinci- 
dence window r we must multiply (1361) by the unit 



1 = dt6(t - L/v T y)dt'S(t' - L/v T z + t) (38) 
and integrate over dydz. As a result we obtain 

( t j \ dt'T 2 

N c . d N ( k„ Pl )IM = JW (-, ) ^p^. (39) 

After integration over c?t' in the range of the coincidence window r we can put z ~ y, and 
finally get 

iV c « iV ^ p ^exp(-2 2 / 2 )-. (40) 

For comparison of TOF spectrum of two and single detectors it is useful to find ratio of (1401) 
to fl34|) . This ratio is 

N r 

W = -^=dJ2p-fl(y), (41) 

^) = ^yexp(-V). (42) 

The function R(y) is shown in Fig. 3. Its integral / dyR(y) is equal to 0.15. So we can tell that 
the ratio is approximately 

W « O.ldfi- (43) 

as is said in section 2. 



where 
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Figure 3: Dependence of R(y) on y = k/yhvT. 

5 Conclusion 

We have shown that the effect of n-n scattering experiment and spectrum of detected neutrons 
in a single detector can be calculated analytically with the standard scattering theory without 
transformation to center of mass system. Analytically calculated factor F = 0.705 is close 
to that F = 0.83, calculated by Monte Carlo method. The difference can be attributed to 
slightly different spectra of neutrons in the collision volume. In Monte Carlo calculations 
spectrum contained Maxwellian part and epithermal tail, while for analytical calculations we 
used only Maxwellian part. We did not calculated background which is related to scattering 
of neutrons on gas molecules, but we claim that it also can be calculated analytically. One 
of the main conclusions of this paper is that coincidence scheme for this type of experiment 
is absolutely impractical, because the effect becomes so low, that the level of the background 
becomes irrelevant. 
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6 The history of submissions and rejections 

I submitted this paper on March 10 to the same journal J. Phys. G: Nucl. Part. Phys., where 
the two papers[Tj, [2] were published. On March 31 I received the electronic mail from editors 
with subject Final decision on your article from J. Phys. G: Nucl. Part. Phys., which 
meant that no more negotiations are supposed. The letter contained the referee report. I am 
not permitted by arXiv policy to present full content of the report, so I give a paraphrase of it. 

6.1 In his report the referee writes 

that the paper is not worth of further consideration by the journal, because the journal has 
already published articles where neutron spectra were calculated by Monte Carlo techniques. 
Analytical calculations by standard technique are not worth to be published. A new result of the 
paper, which shows that coincidence scheme is not profitable was also obtained by participants 
of the project (their result was not yet published or submitted for publication), therefore it is 
necessary to reject the paper. 
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6.2 My remark not sent to the editor 



In my opinion analytical calculation is much superior to Monte Carlo one. If it was not known 
and became known only later, it is worth its own publication. Even if it gives the same result 
as Monte Carlo one, it is very important, because it proves that Monte Carlo calculations was 
made correctly. No one can check Monte Carlo calculations. Every one can check analytical 
calculations. It is great that for analytical calculations one need nothing except standard theory, 
because, if it were necessary to make some new assumption, the merit of analytical calculation 
would become doubtful. 

Referee also pointed out, that the papers [1,2] contain calculation of signal and background, 
while I calculated only signal. To this remark I can say that I am able to calculate the back- 
ground too. However this analytical calculation is more laborious and why to do it, if it had 
already been calculated by the Monte-Carlo method, and analytical calculation are not worth 
of a "stand-alone" publication? 

After getting this reply I reconsidered my article and improved it. But in essence it remained 
the same. So the referee, if he were able to read the improved version, would have no reason 
to change his report. 

6.3 The end of the story 

On April 3 I submitted the paper to Yad.Phys. (Russian Nuclear Physics), and on April 29 
I received the referee report, which approved the paper and contained some comments that 
helped me to improve it even further. So I want to express my gratitude to him. 

6.4 The story continues 

After correction of the article I submitted it again to Yad. Phys. but after few days I obtained 
the letter, where editors asked me to delete section 6. I replied why? Whether it is not my right 
to publish everything, which is related to the problem? More over the referee did not require 
to omit this part. But in the next letter the editors informed me that the referee considered 
the deletion of the section 6 as self-evident. I sent to the editorial board my arguments why to 
accept my paper with the section 6. They are: 

1. It is nasty to forbid something which can be permitted. 

2. It is nasty to apply power where it is useless. 

3. Publication of the section 6 is harmful only for editorial board of the J. Phys. G, and 
their referee, but is very profitable for the whole scientific community, because it shows 
that irresponsible referee reports will be published and it is a real punishment for them. 

4. Any decision of the Yad. Phys. editors will be historical one, because I shall publish 
everything in the ArXiv, but the positive decision will demonstrate that the editorial 
board agree with me and takes the responsibility for future reports of its referees. 

5. The section 6 is not irrelevant to the content of the paper. It rises the important ques- 
tion: whether analytical calculations merit publication as a stand-alone article or not, if 
everything can be calculated numerically by, say, Monte-Carlo method. 

Notwithstanding of my arguments the verdict was - to delete the section 6. One of the vice 
chief editors (he is not anonymous and advised me not to reveal his name) wrote me that he 
CANNOT publish because it contradicts to the LAW of GENRE. I asked him which article 
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of the GENRE LAW does he refers? One of the Journals publishes my articles with referee 
reports. The editorial board does not think that it contradicts the GENRE LAW and I take 
off my hat to their editors. However, since there are no arguments except the claim that the 
Journal CANNOT publish the section 6, I consider it as a demonstration of power to which I 
am to obey. I delete the section, but I shall publish the full article in the ArXiv, and there I 
shall explain why the section 6 is excluded from the article published in Yad.Phys. 

6.5 An attempt to publish JINR preprint 

I wanted also to publish my paper with section 6 as a preprint of JINR. The chief of publishing 
department told me that she has no right to publish such a section without approval by Scientific 
Secretary (SS). I applied to SS for permission, but in vain. Then I appealed to our director (D) 
with the same arguments as above. 

-Stop! — said one of my friends to whom I told the story. - I can predict the end of it. 
The D will sent your article again to SS. The SS will return it to D with a note: "I consider it 
not appropriate" , and the D will write "I agree" and sign such a resolution! 

I was astonished how smart was my friend! 
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